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Abstract 

We know that giant planets played a crucial role in the making of 
our Solar System. The discovery of giant planets orbiting other stars 
is a formidable opportunity to learn more about these objects, what is 
their composition, how various processes influence their structure and 
evolution, and most importantly how they form. Jupiter, Saturn, Uranus 
and Neptune can be studied in detail, mostly from close spacecraft flybys. 
We can infer that they are all enriched in heavy elements compared to 
the Sun, with the relative global enrichments increasing with distance 
to the Sun. We can also infer that they possess dense cores of varied 
masses. The intercomparison of presently caracterised extrasolar giant 
planets show that they are also mainly made of hydrogen and helium, but 
that they either have significantly different amounts of heavy elements, 
or have had different orbital evolutions, or both. Hence, many questions 
remain and are to be answered for significant progresses on the origins of 
planets. 
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— Pourquoi Vazur muet et I'espace insondable? 
Pourquoi les astres d'or fourmillant comme un sable? 

Arthur Rimbaud — Soleil et chair 

1 Introduction 

Looking at a starry sky, it is quite vertiginous to think that we are at one of 
these special epochs in history: Just before, we only knew of the planets in our 
Solar System. Now, more than 150 giant planets are known to orbit solar-like 
stars. Our giant planets, Jupiter, Saturn, Uranus and Neptune are no longer 
the only ones that we can characterize. We now know of six extrasolar giant 
planets transiting in front of their stars, for which we can measure with a fair 
accuracy their mass and radius. We lie on the verge of a true revolution: With 
ground-based and future space-based transit search programs, we should soon 
be able to detect and characterize many tens, probably hundreds of planets 
orbiting their stars, with the hope of inferring their composition and hence the 
mechanisms responsible for the formation of planets. 

It is a daunting task too, because we should expect that, like for the planets 
in our Solar System, a rich variety of giant planets is found, with different 
compositions, different histories and a number of new or unexpected physical 
mecanisms at work. We will have to classify observations, test theories, and be 
aware that although simplicity is appealing, it is not always what Nature has 
in store for us. 

This review aims at providing a synthetic approach to the problems posed 
by "old" and "new" giant planets, in the Solar System and outside. It updates 
a previous review by Stevenson (1982), and expands on the review by Hubbard 
et al. (2002) by focusing on the mass-radius relations and compositions of giant 
planets. In Section 2, we see how to construct interior and evolution models 
of giant planets. Section 3 is devoted to our giant planets, what we can infer 
from observations, and the questions that remain. I then turn to the new field 
of extrasolar giant planets, focusing on the close-in, transiting "Pegasi planets" 
(also called "hot Jupiters" ) . The last section is an attempt to summarize some of 
the known facts concerning giant planets and provide a few expected milestones 
for future studies. 

2 The calculation of interior models 
2.1 A simple model 

To tell our story, I will use a simple model, based on the following assumptions: 

1. Giant planets are made of a fluid envelope and possibly a dense central 
core of about ~ 15 M s (Earth masses); 

2. The envelope is mostly made of hydrogen and helium and trace species 
(heavy elements); The core is made of an unknown combination of refrac- 
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tory material ( "rocks" ) and more volatile species ( "ices" including molec- 
ular species such as H20 water, CH4 methane and NH3 ammonia in the 
fluid state); 

3. Contrary to solid planets, viscosity is negligible throughout; 

4. In most cases, rotation and magnetic fields can be neglected; 

5. Giant planets were formed from an extended, high-entropy, high-luminosity 
state. 

These assumptions can only be justified a posteriori: they are the result of our 
knowledge of observed giant planets and of inferences about the mechanisms that 
led to their formation. We will see how this simple model predicts the global 
properties of giant planets between 1/20 to 20 Mj (about 15 and 6000 M e ), and 
how they compare with observations. 



2.2 Basic equations 

As a consequence of our assumptions, the structure and evolution of a giant 
planet is governed by the following hydrostatic, thermodynamic, mass conser- 
vation and energy conservation equations: 
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where P is the pressure, p the density, and g = Gm/r 2 the gravity (m is the 
mass, r the radius and G the gravitational constant). The temperature gradient 
Vt = (d In T/d In P) depends on the process by which the internal heat is 
transported. L is the intrinsic luminosity, t the time, S the specific entropy (per 
unit mass), and e accounts for the sources of energy due e.g. to radioactivity 
or more importantly nuclear reactions. Generally it is a good approximation to 
assume e <~ for objects less massive than ~ 13 Mj, i.e. too cold to even burn 
deuterium (but we will see that in certain conditions this term may be useful, 
even for low mass planets). 

The boundary condition at the center is trivial: r = 0; (m = 0, L = 0). 
The external boundary condition is more difficult to obtain because it depends 
on how energy is transported in the atmosphere. One possibility is to use the 
Eddington approximation, and to write (e.g. Chandrasckhar 1939): r = R; 
(T = T e g, P = 2/3. g/«), where T e g is the effective temperature, and k is the 
opacity in cm 2 g _1 . Note for example that in the case of Jupiter T e s = 124 K, 
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g = 2600 cm s" 2 and « « 5 x 10~ 2 (P/1 bar) cm 2 g" 1 . This implies F ~ 0.2 bar, 
which is actually close to Jupiter's tropopause, where T « 110 K. 

More generally, one has to use an atmospheric model relating the tempera- 
ture and pressure at a given level to the radius R, intrinsic luminosity L and in- 
coming stellar luminosity L* p : r = R; (T Q = T (R, L, L* p ), Pq = Po(R, L, L* p )). 
P Q is chosen to satisfy the condition that the corresponding optical depth at that 
level should be much larger than unity. If the stellar flux is absorbed mostly in 
a convective zone, then the problem can be simplified by using Tq(R, L, L* p ) rj 
T (R, L + L^p, 0) (e.g. Hubbard 1977). An example of such a model is described 
by Saumon et al. (1996) and Hubbard et al. (2002) and is used hereafter to 
model the planets in the low irradiation limit. 

2.3 High pressure physics & equations of state 

In terms of pressures and temperatures, the interiors of giant planets lie in a 
region for which accurate equations of state (EOS) are extremely difficult to 
calculate. This is because both molecules, atoms and ions can coexist, in a fluid 
that is partially degenerate (free electrons have energies that are determined 
both by quantic and thermal effects) and partially coupled (coulombian inter- 
actions between ions are not dominant but must be taken into account). The 
presence of many elements and their possible interactions further complicate 
matters. For lack of space, this section will mostly focus on hydrogen whose 
EOS has seen the most important developments in recent years. A phase di- 
agram of hydrogen (fig. ^| illustrates some of the important phenomena that 
occur in giant planets. 

The photospheres of giant planets are generally relatively cold (50 to 3000 K) 
and at low pressure (0.1 to lObar, or 10 4 to 10 6 Pa), so that hydrogen is in 
molecular form and the perfect gas conditions apply. As one goes deeper into 
the interior hydrogen and helium progressively become fluid. (The perfect gas 
relation tends to underestimate the pressure by 10% or more when the density 
becomes larger than about 0.02 gem -3 (P > lkbar in the case of Jupiter)). 

Characteristic interior pressures are considerably larger however: as implied 
by Eqs. HI and P c ~ GM 2 /R i , of the order of 10-100 Mbar for Jupiter and 
Saturn. At these pressures and the corresponding densities, the Fermi tem- 
perature Tp is larger than 10 5 K. This implies that electrons are degenerate. 
Figure^shows that inside Jupiter, Saturn, HD209458b, but also for giant plan- 
ets in general for most of their history, the degeneracy parameter 6 = T/Tp 
is between 0.1 and 0.03. Therefore, the energy of electrons in the interior is 
expected to be only slightly larger than their non-relativistic, fully degenerate 
limit: u c > 3/5/cTf = 15.6 (p/fic) 2 ^ 3 eV, where k is Boltzmann's constant, 
/i is the number of electrons per nucleon and p is the density in gcm~ 3 . For 
pure hydrogen, when the density reaches ~ 0.8 gem -3 , the average energy of 
electrons becomes larger than hydrogen's ionization potential, even at zero tem- 
perature: hydrogen pressure-ionizes and becomes metallic. This molecular to 
metallic transition occurs near Mbar pressures, but exactly how this happens 
remains unclear because of the complex interplay of thermal, coulombian, and 
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Figure 1: Phase diagram for hydrogen with the main phase transitions occurring 
in the fluid or gas phase. The temperature-pressure profiles for Jupiter, Saturn, 
Uranus, Neptune, and HD209458 b are shown. The dashed nearly vertical 
line near 1 Mbar is indicative of the molecular to metallic transition (here it 
represents the so-called plasma phase transition as calculated by Saumon et al. 
2000). The region in which hydrogen is in solid phase (Datchi et al. 2000; 
Gregoryanz et al. 2003) is represented as a dashed area. The three phases 
(I, II, III) of solid hydrogen are shown (see Mao & Hemley 1994). Values of the 
degeneracy parameter 8 are indicated as dotted lines to the upper right corner 
of the figure. 
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degeneracy effects (in particular, whether hydrogen metallizes into an atomic 
state H + — as suggested in Fig. 1 — or first metallizes in the molecular state 
H2 remains to be clarified). 

Recent laboratory measurements on fluid deuterium have been able to reach 
pressures above <; 1 Mbar, and provide new data in a region where the EOS 
remains most uncertain. Gas-guns experiments have been able to measure the 
reshock temperature (Holmes et al. 1995), near T ~ 5000 K, P ~ 0.8 Mbar, 
and a rise in the conductivity of molecular hydrogen up to T ~ 3000 K, P ~ 
1.4 Mbar, a sign that metallicity may have been reached (Weir et al. 1996). The 
following few years have seen the development of laser-induced shock compres- 
sion (Da Silva et al. 1997, Collins et al. 1998), pulsed-power shock compression 
(Knudson et al. 2002, 2004), and convergent shock wave experiments (Belov 
et al. 2002; Boriskov et al. 2003) in a high-pressure (P = 0.3 - 4 Mbar) high- 
temperature (T ~ 6000 — 10 5 K) regime. Unfortunately, experimental results 
along the principal Hugoniot of deuterium do not agree in this pressure range. 
Laser compression data give a maximum compression of ~ 6 while both the 
pulsed-power compression experiments and the convergent shock wave exper- 
iments find a value of ~ 4 . Models that are partly calibrated with experi- 
mental data (Saumon, Chabrier & Van Horn 1995; Ross 1998; Saumon et al. 
2000, Ross & Yang 2001) obtain a generally good agreement with the laser- 
compression data. However, the fact that independant models based on first 
principles (Militzer & Ceperley 2001; Desjarlais 2003; Bonev et al. 2004) yield 
low compressions strongly favors this solution. 

The question of the existence of a first-order molecular to metallic transition 
of hydrogen (i.e. both molecular dissociation and ionisation occur simultane- 
ously and discontinuously at the so-called plasma phase transition, or PPT) 
remains however. The critical line shown in fig. ^ corresponds to calculations 
by Saumon et al. (2000), but may be caused by artefacts in the free energy 
calculation. Recent Density Functional Theory (DFT) simulations by Bonev 
et al. (2004) indicate the possibility of a first order liquid-liquid transition but 
other path-integral calculations (Militzer & Ceperley 2001) do not. It is crucial 
to assess the existence of such a PPT because it would affect both convection 
and chemical composition in the giant planets. 

A clear result from fig. Q]at least is that, as first shown by Hubbard (1968), 
the interiors of the hydrogen-helium giant planets are fluid, whatever their 
age (an isolated Jupiter should begin partial solidification only after at least 
~ 10 3 Ga of evolution). For Uranus and Neptune, the situation is actually more 
complex because at large pressures they are not expected to contain hydrogen, 
but numerical simulations show that ices in their interior should be fluid as well 
(Cavazzoni et al. 1999). 

Models of the interiors of giant planets require thcrmodynamically consistent 
EOSs calculated over the entire domain of pressure and temperature spanned by 
the planets during their evolution. Elements other than hydrogen, most impor- 
tantly helium, should be consistently included. Such a calculation is a daunting 
task, and the only recent attempt at such an astrophysical EOS for substellar 
objects is that by Saumon et al. (1995). Another set of EOSs reproducing cither 
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the high- or low-compression results was calculated by Saumon & Guillot (2004) 
specifically for the calculation of present-day models of Jupiter and Saturn. 

These EOSs have so far included other elements (including helium), only in 
a very approximative way, i.e. with EOSs for helium and heavy elements that 
are based on interpolations between somewhat ideal regimes, using an additive 
volume law, and neglecting the possibility of existence of phase separations (see 
Hubbard et al. 2002 and Guillot et al. 2004 for further discussions). 

2.4 Heat transport 

Giant planets possess hot interiors, implying that a relatively large amount of 
energy has to be transported from the deep regions of the planets to their surface. 
This can either be done by radiation, conduction, or, if these processes are not 
sufficient, by convection. Convection is generally ensured by the rapid rise of the 
opacity with increasing pressure and temperature. At pressures of a bar or more 
and relatively low temperatures (less than 1000 K), the three dominant sources 
of opacities are water, methane and collision-induced absorption by hydrogen 
molecules. 

However, in the intermediate temperature range between ~ 1200 and 1500 K, 
the Rosscland opacity due to the hydrogen and helium absorption behaves dif- 
ferently: the absorption at any given wavelength increases with density, but 
because the temperature also rises, the photons are emitted at shorter wave- 
lengths, where the monochromatic absorption is smaller. As a consequence, the 
opacity can decrease. This was shown by Guillot et al. (1994) to potentially 
lead to the presence of a deep radiative zone in the interiors of Jupiter, Saturn 
and Uranus. 

This problem must however be reanalyzed in the light of recent observations 
and analyses of brown dwarfs. Their spectra show unexpectedly wide sodium 
and potassium absorption lines (see Burrows, Mar ley & Sharp 2000), in spectral 
regions where hydrogen, helium, water, methane and ammonia are relatively 
transparent. It thus appears that the added contribution of these elements (if 
they are indeed present) would wipe out any radiative region at these levels 
(Guillot et al. 2004). 

At temperatures above 1500 ~ 2000 K two important sources of opacity 
appear: (i) the rising number of electrons greatly enhances the absorption of 
H^~ and H~; (ii) TiO, a very strong absorber at visible wavelengths is freed 
by the vaporization of CaTiOa. Again, the opacity rises rapidly which ensures 
a convective transport of the heat. Still deeper, conduction by free electrons 
becomes more efficient, but the densities are found not to be high enough for 
this process to be significant, except perhaps near the central core (see Hubbard 
1968; Stevenson & Salpeter 1977). 

However, because irradiated giant planets do develop a radiative zone, Rossc- 
land opacity tables covering the proper range of temperatures and pressures are 
needed. A pure hydrogen-helium mixture table has been calculated by Lenzuni 
et al. (1991). Opacities for solar composition including grains are available 
from Alexander & Ferguson (1994), but they do not include alkali metals and 
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up-to-date data on water, methane and TiO absorption. Guillot (1999a) pro- 
vides a grain-free, alkali-free table which is limited to low-temperature regimes. 
The calculations hereafter use opacities provided by F. Allard on the basis of 
calculations for brown dwarfs of solar composition, including grains and alkali 
metals (Allard et al. 2001). 



2.5 The contraction and cooling histories of giant planets 

The interiors of giant planets is expected to evolve with time from a high entropy, 
high 9 value, hot initial state to a low entropy, low 8, cold degenerate state. The 
essential physics behind can be derived from the well-known virial theorem and 
the energy conservation which link the planet's internal energy E;, gravitational 
energy E g and luminosity through: 

m + E e = 0, (5) 

* - 

where £ = J Q M 3(P/ p)dmj J Q M udra «< 3P/ pu > and u is the specific internal 
energy. For a diatomic perfect gas, £ = 3.2; for fully-degenerate non-relativistic 
electrons, £ = 2. 

Thus, for a giant planet or brown dwarf beginning its life mostly as a perfect 
H2 gas, two third of the energy gained by contraction is radiated away, one 
third being used to increase E;. The internal energy being proportional to the 
temperature, the effect is to heat up the planet. This represents the slightly 
counter-intuitive but well known effect that a star or giant planet initially heats 
up while radiating a significant luminosity. 

Let us now move further in the evolution, when the contraction has pro- 
ceeded to a point where the electrons have become degenerate. For simplicity, 
I will ignore Coulombian interactions and exchange terms, and assume that 
the internal energy can be written as E\ = E c \ + Ei on , and that furthermore 
i?ci 3> -E-ion {0 is small). Because £ « 2, we know that half of the gravitational 
potential energy is radiated away and half of it goes into internal energy. The 
problem is to decide how this energy is split into an electronic and an ionic part. 
The gravitational energy changes with some average value of the interior den- 
sity as £g oc 1/R oc p 1 ^ 3 . The energy of the degenerate electrons is essentially 
the Fermi energy: E c \ oc p 2 / 3 . Therefore, E c \ w 2(E e /E g )E g . Using the virial 
theorem, this yields: 

E « —Eg w 2L (7) 
L w -E ion cx-T. (8) 

The gravitational energy lost is entirely absorbed by the degenerate electrons, 
and the observed luminosity is due to the thermal cooling of the ions. 

Several simplifications limit the applicability of this result (that would be 
valid in the white dwarf regime). In particular, the coulombian and exchange 
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terms in the EOS introduce negative contributions that cannot be neglected. 
However, the approach is useful to grasp how the evolution proceeds: in its 
very early stages, the planet is very compressible. It follows a standard Kelvin- 
Helmoltz contraction. When degeneracy sets in, the compressibility becomes 
much smaller (qT ~ 0.1, where a is the coefficient of thermal expansion), and 
the planet gets its luminosity mostly from the thermal cooling of the ions. The 
luminosity can be written in terms of a modified Kelvin- Helmoltz formula: 

GM 2 

L-r,—, (9) 

where r is the age, and 77 is a factor that hides most of the complex physics. 
In the approximation that Coulombian and exchange terms can be neglected, 
rj w 0/(9+1). The poor compressibility of giant planets in their mature evolution 
stages imply that 77 <C 1: the luminosity is not obtained from the entire gravi- 
tational potential, but from the much more limited reservoir constituted by the 
thermal internal energy. Equation [5] shows that to first order, logL oc — logr: 
very little time is spent at high luminosity values. In other words, the problem 
is (in most cases) weakly sensitive to initial conditions. 




Figure 2: Luminosity versus mass for giant planets after 4.5 Ga of evolution 
compared to measured values for our four giant planets (including the significant 
uncertainty on Uranus' luminosity). The lines correspond to: H+He: a pure 
hydrogen-helium composition with a helium mass mixing ratio Y — 0.25; (a): a 
model with Y = 0.30 and a 15M® core; (b): the same model but with Y = 0.36. 
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Figure[21shows calculated luminosities in the framework of our simple model. 
Compared to J5J, calculated luminosities are consistent with 77 ks 0.01 to 0.03. 
The lower luminosities obtained in the presence of a core and of more heavy 
elements are due to an earlier contraction, and quicker loss of the internal heat. 
As model (b) would be appropropriate to explain Saturn's radius (see next 
section), it can be seen that the planet emits more heat than predicted by 
homogeneous contraction models. The cases of Uranus and Neptune is more 
complex and cannot be directly compared with the models in fig.[21which neglect 
the thermal heat content of the central core. 



2.6 Mass-radius relation 




Figure 3: Radius versus mass for giant planets after 4.5 Ga of evolution com- 
pared to measured values for our four giant planets and four known extrasolar 
planets. As in fig. [3 the lines correspond to: H+Hc: a pure, Y — 0.25, hydrogen- 
helium composition (Y=0.25); (a): a model with Y = 0.30 and a 15 core; 
(b): the same model but with Y = 0.36. An approximate mass-radius relation 
for zero-temperature water and olivine planets is shown as dashed and dash- 
dotted lines, respectively (Courtesy of W.B. Hubbard). 

The relation between mass and radius has very fundamental astrophysical 
applications. Most importantly it allows one to infer the gross composition of an 
object from a measurement of its mass and radius. This is especially relevant in 
the context of the discovery of extrasolar planets with both radial velocimetry 
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and the transit method, as the two techniques yield relatively accurate deter- 
mination of M and R. 

Figure|3shows the mass-radius relation for isolated or nearly-isolated gaseous 
planets, based on our simple simple model and various assumption on their com- 
position. The curves have a local maximum near 4Mj: at small masses, the 
compression is rather small so that the radius increases with mass. At large 
masses, degeneracy sets in and the radius decreases with mass. 

This can be understood on the basis of polytropic models based on the 
assumption that P = Kp 1+1 / n , where K and n are constants. Because of 
degeneracy, a planet of large mass will tend to have n — > 1.5, while a planet a 
smaller mass will be less compressible (n — ► 0). Indeed, it can be shown that in 
their inner 70 to 80% in radius isolated planets of 10, 1 and 0.1 Mj have n = 1.3, 
1.0 and 0.6, respectively. From polytropic equations (e.g. Chandrasekhar 1939): 

RocK^M^. (10) 

Assuming that K is independant of mass, one gets R oc M 016 , M°, and M -018 
for M = 10, 1 and 0.1 Mj, respectively, in relatively good agreement with fig. 01 
(the small discrepancies are due to the fact that the intrinsic luminosity and 
hence K depend on the mass considered). 

Figure [3] shows already that the planets in our Solar System are not made 
of pure hydrogen and helium: their radii lie below that predicted for Y = 0.25 
objects. Indeed, Jupiter, Saturn, and the two ice-giants Uranus and Neptune 
contain a growing proportion of heavy elements. The theoretical curves for 
olivine and ice planets predict even smaller radii however: even Uranus and 
Neptune contain 10 to 20% of their mass as hydrogen and helium. 

The extrasolar planets detected so far (see table OH hereafter) all lie above 
the pure hydrogen-helium curve. This is due to the fact that these planets 
have their evolutions dominated by the intense stellar irradiation they receive. 
Thermal effects are no longer negligible: Using the Eddington approximation, 
assuming k oc P and a perfect gas relation in the atmosphere, one can show 
that K oc (M/R 2 )- 1 / 2n and that therefore R oc . With n = 1, one finds 

R oc M -1 / 2 . Strongly irradiated hydrogen-helium planets of small masses are 
hence expected to have the largest radii which qualitatively explain the positions 
of the extrasolar planets in fig. Note that this estimate implicitly assumes 
that n is constant throughout the planet. The real situation is more complex 
because of the growth of a deep radiative region in most irradiated planets, and 
because of structural changes between the degenerate interior and the perfect 
gas atmosphere. 

2.7 Rotation and the figures of planets 

The mass and radius of a planet informs us on its global composition. Because 
planets are also rotating, one is allowed to obtain more information on their deep 
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interior structure. The hydrostatic equation becomes more complex however: 



^-v(o///*l [ A')-nx(ox,), (ii) 

where fi is the rotation vector. The resolution of eq. Ijlll) is a complex problem. 
It can however be somewhat simplified by assuming that |fi| = lu is such that the 
centrifugal force can be derived from a potential. The hydrostatic equilibrium 
then writes VP = pS7U, and the figure of the rotating planet is then defined by 
the U = cte level surface. 

One can show (e.g. Zharkov & Trubitsyn 1978) that the hydrostatic equation 
of a fluid planet can then be written in terms of the mean radius f (the radius 
of a sphere containing the same volume as that enclosed by the considered 
equipotential surface) : 

ldP Gm 2 2 GM 

- 7F = -— + -w 2 r+^ -F<p u , 12 
p or r z 3 R A 

where M and R are the total mass and mean radius of the planet, and <p u is 
a slowly varying function of f. (In the case of Jupiter, ip^ varies from about 
2 x 10 -3 at the center to 4 x 10~ 3 at the surface.) Equations lj 214(1 remain the 
same with the hypothesis that the level surfaces for the pressure, temperature, 
and luminosity are equipotentials. The significance of rotation is measured by 
the ratio of the centrifugal acceleration to the gravity: 

uj 2 RL 

« = w- < 13 > 

The external gravitational potential of the planet is (assuming hydrostatic 
equilibrium) : 



V exb (r,cos6) = 



(14) 



where the coefficients J 2n are the planet's gravitational moments, and the P2„ 
are Legendre polynomials. The J's can be measured by a spacecraft coming 
close to the planet, preferably on a polar orbit. Together with the mass, this 
provides a constraint on the interior density profile (see Zharkov & Trubitsyn 
1974): 

M = [[[ p(r,8)d 3 r, 



hi = -jjjp- P (r,e)r 2i P 2i (cos 6)d 3 r, 

where dr is a volume element and the integrals are performed over the entire 
volume of the planet. 
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Figure|3|shows how the different layers inside a planet contribute to the mass 
and the gravitational moments. The figure applies to Jupiter, but would remain 
very similar for other planets. Measured gravitational moments thus provide 
information on the external levels of a planet. It is only indirectly, through 
the constraints on the outer envelope that the presence of a central core can 
be infered. As a consequence, it is impossible to determine this core's state 
(liquid or solid), structure (differentiated, partially mixed with the envelope) 
and composition (rock, ice, helium...). 




R/Rtot 



Figure 4: Contribution of the level radii to the gravitational moments of Jupiter. 
Jo is equivalent to the planet's mass. The small discontinuities are caused by the 
following transitions, from left to right: core/envelope, helium rich/helium poor 
(metallic/molecular). Diamonds indicate the median radius for each moment. 

For planets outside the solar system, although measuring their gravitational 
potential is utopic, their oblateness may be reachable with future space transit 
observations (Seager & Hui 2002). Since the oblateness e is, to first order, 
proportionnal to q: 

R cq I 3 . I 



oA 2 + - q (15) 



i? cq — R po \ \ 2 2 

(where A 2 = J2I1 ~ 0.1 to 0.2), it may be possible to obtain their rotation rate, 
or with a rotation measured from another method, a first constraint on their 
interior structure. 
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3 Jupiter, Saturn, Uranus and Neptune 



3.1 Main observational data 

The mass of the giant planets can be obtained with great accuracy from the 
observation of the motions of their natural satellites: 317.834, 95.161, 14.538 
and 17.148 times the mass of the Earth (1M 9 = 5.97369 x 10 27 g) for Jupiter, 
Saturn, Uranus and Neptune, respectively. The more precise determination of 
their gravity fields listed in table have been obtained by the Pioneer and 
Voyager space missions. 



Table 1: Characteristics of the gravity fields and radii 





Jupiter 


Saturn 


Uranus 


Neptune 


M x 10 _2y [g] 


18.986112(15) a 


5.684640(30)° 


0.8683205(34)° 


1.0243542(31) d 


i? cq x 1(T 9 [cm] 


7.1492(4)° 


6.0268(4) f 


2.5559(4) g 


2.4766(15) g 


i?poi x 1CT 9 [cm] 


6.6854(10)° 
6.9894(6) h 


5.4364(10) f 
5.8210(6) h 


2.4973(20) g 


2.4342(30) g 


R x 1CT 9 [cm] 


2.5364(10)' 


2.4625(20)' 


P [gem -3 ] 


1.3275(4) 


0.6880(2) 


1.2704(15) 


1.6377(40) 
0.3539(10) d 
-0.28(22) d 


J 2 x 10 2 


1.4697(l) a 


1.6332(10) b 
-9.19(40) b 
1.04(50)° 


0.35160(32)° 


Ja x 10 4 


-5.84(5) a 


-0.354(41)° 


Ja x 10 4 


0.31(20) a 






P u x 1(T 4 [s] 


3.57297(41) j 


3.83577(47) j 


6.206(4) k 


5.800(20)' 


Q 


0.08923(5) 


0.15491(10) 


0.02951(5) 


0.02609(23) 


C/MR 2 C<1 


0.258 


0.220 


0.230 


0.241 



The numbers in parentheses are the uncertainty in the last digits of the given 

value. The value of the gravitational constant used to calculate the masses of 

Jupiter and Saturn is G = 6.67259 x 10~ 8 dyn.cm^g" 1 (Cohen & Taylor, 1987). 

a Campbell & Synott (1985) 

b Campbell & Anderson (1989) 

° Anderson et al. (1987) 

d Tyler et al. (1989) 

° Lindal et al. (1981) 

f Lindal et al. (1985) 

g Lindal (1992) 

From 4th order figure theory 
1 (27? cq +-R P oi)/3 (Clairaut's approximation) 
j Davies et al. (1986) 
k Warwick et al. (1986) 
1 Warwick et al. (1989) 

Table^also indicates the radii obtained with the greatest accuracy by radio- 
occultation experiments. By convention, these radii and gravitational moments 
correspond to the 1 bar pressure level. The rotation periods are measured from 
the variations of the planets' magnetic fields (system III) and are believed to 
be tied to the interior rotation. The giant planets are relatively fast rotators, 
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with periods of about 10 hours for Jupiter and Saturn, and about 17 hours for 
Uranus and Neptune. The fact that this fast rotation visibly affects the figure 
(shape) of these planets is seen by the significant difference between the polar 
and equatorial radii. 

A first result obtained from the masses and radii (using the planets' mean 
radii, as defined in section l2~7|l indicated in Tablets the fact that these planets 
have low densities. These densities are similar, but considering that compres- 
sion strongly increases with mass, one is led to a sub-classification between the 
hydrogen-helium giant planets Jupiter and Saturn, and the "ice giants" Uranus 
and Neptune. 

The values of the axial moment of inertia C have been calculated using 
the Radau-Darwin approximation (Zharkov k, Trubitsyn 1978). Our four giant 
planets all have an axial moment of inertia substantially lower than the value 
for a sphere of uniform density, i.e. 2/5 MR 2 , indicating that they have dense 
central regions. This does not necessarily mean that they possess a core, but 
simply that the density profile departs significantly from a uniform value. 



Table 2: Energy balance as determined from Voyager IRIS data a . 





Jupiter 


Saturn 


Uranus 


Neptune 


Absorbed power [10^ erg.s -1 ] 


50.14(248) 


11.14(50) 


0.526(37) 


0.204(19) 


Emitted power [10 23 erg.s -1 ] 


83.65(84) 


19.77(32) 


0.560(11) 


0.534(29) 


Intrinsic power [10 23 erg.s -1 ] 


33.5(26) 


8.63(60) 


0.034(38) 


0.330(35) 


Intrinsic flux [erg.s -1 . cm -2 ] 


5440.(430) 


2010.(140) 


42.(47) 


433.(46) 


Bond albedo [] 


0.343(32) 


0.342(30) 


0.300(49) 


0.290(67) 


Effective temperature [K] 


124.4(3) 


95.0(4) 


59.1(3) 


59.3(8) 


1-bar temperature 1 " [K] 


165.(5) 


135.(5) 


76.(2) 


72.(2) 



After Pearl & Conrath (1991) 



b Lindal (1992) 

Jupiter, Saturn and Neptune are observed to emit significantly more energy 
than they receive from the Sun (see Table 0). The case of Uranus is less clear. 
Its intrinsic heat flux F{ nt is significantly smaller than that of the other giant 
planets. Detailed modeling of its atmosphere however indicate that F- Ult > 
60 erg cm -2 s -1 (Marley & McKay 1999). With this caveat, all four giant planets 
can be said to emit more energy than they receive from the Sun. Hubbard 
(1968) showed in the case of Jupiter that this can be explained simply by the 
progressive contraction and cooling of the planets. 

It should be noted that the 1 bar temperatures listed in table El are retrieved 
from radio-occultation measurements using a helium to hydrogen ratio which, 
at least in the case of Jupiter and Saturn, was shown to be incorrect. The new 
values of Y are found to lead to increased temperatures by ^ 5 K in Jupiter and 
~ 10 K in Saturn (see Guillot 1999a). However, the Galileo probe found a 1 bar 
temperature of 166 K (Seiff et al. 1998), and generally a good agreement with 
the Voyager radio-occultation profile with the wrong He/H2 value. 
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3.2 Atmospheric composition 

The most important components of the atmospheres of our giant planets are 
also among the most difficult to detect: H 2 and He have a zero dipolar moment. 
Also their rotational lines are either weak or broad. On the other hand, lines due 
to electronic transitions correspond to very high altitudes in the atmosphere, 
and bear little information on the structure of the deeper levels. The only 
robust result concerning the abundance of helium in a giant planet is by in situ 
measurement by the Galileo probe in the atmosphere of Jupiter (von Zahn et 
al. 1998). The helium mole fraction (i.e. number of helium atoms over the total 
number of species in a given volume) is qh c = 0.1359 ±0.0027. The helium mass 
mixing ratio Y (i.e. mass of helium atoms over total mass) is constrained by its 
ratio over hydrogen, X: Y/(X + Y) = 0.238 ± 0.05. This ratio is by coincidence 
that found in the Sun's atmosphere, but because of helium sedimentation in the 
Sun's radiative zone, it was larger in the protosolar nebula: Y proto — 0.275 ±0.01 
and (X + y) P roto ~ 0.98. Less helium is therefore found in the atmosphere of 
Jupiter than inferred to be present when the planet formed. 

Helium is also found to be depleted compared to the protosolar value in 
Saturn's atmosphere. However, in this case the analysis is complicated by the 
fact that Voyager radio occultations apparently led to a wrong value. The 
current adopted value is now Y — 0.18 — 0.25 (Conrath & Gautier 2000), in 
agreement with values predicted by interior and evolution models (Guillot 1999a; 
Hubbard et al. 1999). Finally, Uranus and Neptune are found to have near- 
protosolar helium mixing ratios, but with considerable uncertainty (Gautier & 
Owen 1989). 

The abundance of "heavy elements" , i.e. elements other than hydrogen and 
helium, bears crucial information for the understanding of the processes that 
led to the formation of these planets. Again, the most precise measurements 
are for Jupiter, thanks to the Galileo probe. As shown by fig. |SJ most of the 
heavy elements are enriched by a factor 2 to 4 compared to the solar abundance 
(Niemann et al. 1998; Owen et al. 1999). One exception is neon, but an explana- 
tion is its capture by the falling helium droplets (Roustlon & Stevenson 1995). 
Another exception is water, but this molecule is affected by meteorological pro- 
cesses, and the probe was shown to have fallen into a dry region of Jupiter's 
atmosphere. There are strong indications that its abundance is at least solar. 
Possible very high interior abundances (~ 10 times the solar value) have also 
been suggested as a scenario to explain the delivery of heavy elements to the 
planet (Gautier et al. 2001, Hersant et al. 2004). 

In the case of Saturn, both carbon in the form of methane and nitrogen as 
ammonia appear to be significantly enriched, but with large error bars (Atreya 
et al. 2003). In Uranus and Neptune, methane is probably between 30 and 60 
times the solar value (Gautier & Owen 1989; Hersant et al. 2004). 
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Figure 5: Elemental abundances measured in the tropospheres of Jupiter (cir- 
cles) and Saturn (squares) in units of their abundances in the protosolar nebula. 
The elemental abundances for Jupiter are derived from the in situ measurements 
of the Galileo probe (e.g. Mahaffy et al. 2000; Atreya et al. 2003). Note that the 
oxygen abundance is considered to be a minimum value due to meteorological 
effects (Roos-Serote et al. 2004). The abundances for Saturn are spectroscopic 
determination (Atreya et al. 2003 and references therein). The solar or pro- 
tosolar abundances used as a reference are from Lodders (2003). The arrows 
show how abundances are affected by changing the reference protosolar abun- 
dances from those of Anders & Grevesse (1989) to those of Lodders (2003). 
The horizontal dotted lines indicate the locus of a uniform 2- and 4-times solar 
enrichment in all elements except helium and neon, respectively. 
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Ices + Rocks 
core ? 

Jupiter Saturn 

Figure 6: Schematic representation of the interiors of Jupiter and Saturn. The 
range of temperatures is estimated using homogeneous models and including a 
possible radiative zone indicated by the hashed regions. Helium mass mixing 
ratios Y are indicated. The size of the central rock and ice cores of Jupiter and 
Saturn is very uncertain (see text). In the case of Saturn, the inhomogeneous 
region may extend down all the way to the core which would imply the formation 
of a helium core. [Adapted from Guillot 1999b]. 
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3.3 Interior models: Jupiter and Saturn 

As illustrated by fig. El the simplest interior models of Jupiter and Saturn match- 
ing all observational constraints assume the presence of three main layers: (i) 
an outer hydrogen-helium envelope, whose global composition is that of the 
deep atmosphere; (ii) an inner hydrogen-helium envelope, enriched in helium 
because the whole planet has to fit the H/He protosolar value; (iii) a central 
dense core. Because the planets are believed to be mostly convective, these 
regions are expected to be globally homogeneous. (Many interesting thermo- 
chemical transformations take place in the deep atmosphere, but they are of 
little concern to us). 

A large part of the uncertainty in the models lies in the existence and location 
of an inhomogeneous region in which helium separates from hydrogen to form 
helium-rich droplets that fall deeper into the planet due to their larger density. 
Models have generally assumed this region to be relatively narrow, because 
helium was thought to be most insoluble in low-pressure metallic hydrogen (e.g. 
Stevenson 1982). However, DFT calculations have indicated that the critical 
temperature for helium demixing may rise with pressure (Pfaffcnzeller et al. 
1995), presumably in the regime where hydrogen is only partially ionized and 
bound states remain. This opens up the possibility that the inhomogeneous 
regions may be more extended. In particular, in the case of Saturn, Fortney 
& Hubbard (2003) have shown that explaining Saturn's age may require that 
helium fall all the way to the core, thereby yielding the formation of a helium 
core (or of a helium shell around a rock or ice core). 

With these caveats, the three-layer models can be used as a useful guidance 
to a necessarily hypothetical ensemble of allowed structures and compositions 
of Jupiter and Saturn. Figure shows such an ensemble for Jupiter, based on 
calculations by Saumon & Guillot (2004). The calculations assume that only 
helium is inhomogeneous in the envelope (the abundance of heavy elements is 
supposed to be uniform accross the molecular/metallic hydrogen transition). 
Many sources of uncertainties are included however; among them, the most 
significant are on the equations of state of hydrogen and helium, the uncertain 
values of J4 and Jq, the presence of differential rotation deep inside the planet, 
the location of the helium-poor to helium-rich region, and the uncertain helium 
to hydrogen protosolar ratio. 

These results show that Jupiter's core is smaller than ~ lOM^, and that 
its global composition is pretty much unknown (between 10 to 42 M ffi of heavy 
elements in total) . The models indicate that Jupiter is enriched compared to the 
solar value, particularly with the new, low value of Zq (Lodders 2003) used in 
fig. H This enrichment could be compatible with a global uniform enrichment 
of all species near the atmospheric Galileo values. Alternatively, species like 
oxygen (as mostly water) may be significantly enriched. 

Most of the constraints are derived from the values of the radius (or equiv- 
alently mass) and of J2. The measurement of J4 allows to further narrow the 
ensemble of possible models, and in some cases, to rule out EOS solutions (in 
particular those indicating relatively large core masses, between 10 and 20Mq). 
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M z /Me 

Figure 7: Constraints on Jupiter's interior structure based on Saumon & Guillot 
(2004). The value of the core mass (M core ) is shown in function of the mass of 
heavy elements in the envelope (Mz) for models matching all available observa- 
tional constraints. The dashed region corresponds to models matching the laser 
compression experiments. The plain box corresponds to models matching the 
pulsed p ow6r i md convergent shock compression experiments (see text). Grey 
lines indicate the values of Mz that imply uniform enrichments of the enve- 
lope in heavy elements by factors 2 to 8 times the solar value (Z Q = 0.0149), 
respectively. 
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As discussed in Guillot (1999a) and Saumon & Guillot (2004), most of the uncer- 
tainty in the solution arises because very different hydrogen EOSs are possible. 
The fact that more laboratory and numerical experiments seem to indicate rela- 
tively low-compressions for hydrogen at Mbar pressures points towards smaller 
core masses and a larger amount of heavy elements in the planet (plain box in 
fig. 0). However, this relies on uncertain temperature gradients, because the 
EOSs are based on laboratory data obtained at temperatures higher than those 
relevant to the planetary interiors. 

Results slightly outside the boxes of fig. [7| are possible in the presence of a 
discontinuity of the abundance of heavy elements in the interior. Thus, Guil- 
lot (1999a) found slightly larger core masses (up to 12M®) in the case of the 
Saumon-Chabrier EOS with a first order plasma-phase transition. 




Figure 8: Same as fig. [7|in the case of Saturn. Note that smaller core masses 
could result either from allowing a variation of the abundance of heavy near the 
molecular/metallic transition (Guillot 1999a), or from the presence of a helium 
shell around the core (Fortney & Hubbard 2003). 

In the case of Saturn (fig. El, the solutions depend less on the hydrogen EOS 
because the Mbar pressure region is comparatively smaller. The total amount 
of heavy elements present in the planet can therefore be estimated with a better 
accuracy than for Jupiter. However, because Saturn's metallic region is deeper 
into the planet, it mimics the effect that a central core would have on Ji. If we 
allow for variations in the abundance of heavy elements together with the helium 
discontinuity, then the core mass can become much smaller, and even solutions 
with no core can be found (Guillot 1999a). These solutions depend on the 
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hypothetic phase separation of an abundant species (e.g. water), and generally 
cause an energy problem because of the release of considerable gravitational 
energy. However, another possibility is through the formation of an almost pure 
helium shell around the central core, which could lower the core masses by up 
to 7M,9 (Fortney & Hubbard 2003; Hubbard, personnal communication). 

3.4 Interior models: Uranus and Neptune 




Uranus Neptune 



Figure 9: Schematic representation of the interiors of Uranus and Neptune. 
[Adapted from Guillot 1999b]. 

Although the two planets are relatively similar, fig. already shows that 
Neptune's larger mean density compared to Uranus has to be due to a slightly 
different composition: either more heavy elements compared to hydrogen and 
helium, or a larger rock/ice ratio. The gravitational moments impose that the 
density profiles lie close to that of "ices" (a mixture initially composed of H2O, 
CH4 and NH3, but which rapidly becomes a ionic fluid of uncertain chemical 
composition in the planetary interior), except in the outermost layers, which 
have a density closer to that of hydrogen and helium (Marley et al. 1995; Podolak 
et al. 2000). As illustrated in fig. El three-layer models of Uranus and Neptune 
consisting of a central "rocks" core (magnesium-silicate and iron material), an 
ice layer and a hydrogen-helium gas envelope have been calculated (Podolak et 
al. 1991; Hubbard et al. 1995). 

The fact that models of Uranus assuming homogeneity of each layer and adi- 
abatic temperature profiles fail in reproducing its gravitational moments seem 
to imply that substantial parts of the planetary interior are not homogeneously 
mixed (Podolak et al. 1995). This could explain the fact that Uranus' heat flux 
is so small: its heat would not be allowed to escape to space by convection, but 
through a much slower diffusive process in the regions of high molecular weight 
gradient. Such regions would also be present in Neptune, but much deeper, 
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thus allowing more heat to be transported outward. The existence of these non- 
homogeneous, partially mixed regions arc further confirmed by the fact that if 
hydrogen is supposed to be confined solely to the hydrogen-helium envelope, 
models predict ice/rock ratios of the order of 10 or more, much larger than 
the protosolar value of ~ 2.5. On the other hand, if we impose the constraint 
that the ice/rock ratio is protosolar, the overall composition of both Uranus and 
Neptune is, by mass, about 25% rocks, 60 — 70% ices, and 5 — 15% hydrogen and 
helium (Podolak et al. 1991, 1995; Hubbard et al. 1995). Assuming both ices 
and rocks are present in the envelope, an upper limit to the amount of hydrogen 
and helium present is ~ 4.2 M® for Uranus and ~ 3.2 M ffi for Neptune (Podolak 
et al. 2000). A lower limit of ~ 0.5 M® for both planets can be inferred by 
assuming that hydrogen and helium arc only present in the outer envelope at 
P <, lOOkbar. 

3.5 Are the interiors adiabatic? 

As discussed, the near-adiabaticity of the interiors of the giant planets is a conse- 
quence of the rapid rise of opacities with increasing pressure and temperatures. 
Several exceptions are possible: 

(i) In the "meteorological layer" , the temperature gradient could become 
cither subadiabatic (because of latent heat release and moist convection) or su- 
peradiabatic (because of molecular weight gradients created by condensation 
and precipitation). Locally, a depletion of an efficient radiative absorber (e.g. 
water or methane) could imply that convection is suppressed, either because of 
a lowered radiative gradient, or because sunlight can then be deposited to this 
level. In Uranus and Neptune, a superadiabatic region at P ~ 1 — 2 bar is cor- 
related with methane condensation (Lindal 1992, Guillot 1995). In Jupiter, the 
Galileo probe measured a nearly-adiabatic profile, with a slight static stability 
(N < 0.2 Kkm" 1 ) down to 20 bars (Magalhaes et al. 2002). 

(ii) At the Plasma Phase Transition between molecular and metallic hy- 
drogen, if it exists, with an entropy jump that could be of order 1 fee/baryon 
(Stevenson & Salpeter 1977; Saumon et al. 1995). 

(iii) In the hydrogen-helium phase separation region, where a slow droplet 
formation may inhibit convection and yield a significant superadiabacity (Steven- 
son & Salpeter 1977). 

(iv) Near the core/envelope interface (whether it is abrupt or not) where an 
inhibiting molecular weight gradient occurs and, in the case of Jupiter, conduc- 
tion might play a role. 

(v) Throughout the planets, even though mixing-length arguments predict 
that the superadiabacity is extremely small (<~ 10~ 6 or less), rotation and mag- 
netic fields may increase it, although probably by modest amounts (Stevenson 
1982; See also discussion in Guillot et al. 2004). 
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3.6 What are the ages of our giant planets? 

If we understand something of the formation of our Solar System and of other 
stars, our giant planets should have formed 4.55 Ga ago (e.g. Bodenheimer & 
Lin 2002). The model ages show significant deviations from that value, however. 

In the case of Jupiter, the present radius and luminosity are obtained after 
3.5 to 5.5 Ga of evolution, but most realistic EOSs predict ages above 4.5 Ga 
(Saumon & Guillot 2004). Several processes, among which core erosion, could 
lead to a reduction of that value (Guillot et al. 2004). For Saturn, homogeneous 
evolution models predict ages of order 2 Ga (Stevenson 1982; Saumon et al. 
1992; Guillot et al. 1995). In both planets, the presence of a phase separation 
of helium is likely and would tend to lengthen the cooling. 

The case of Uranus and Neptune is less clear-cut because of the uncertainties 
both on the properties of their atmospheres (in particular their evolution with 
time), and on the global specific heat of material inside. It appears however that 
both planets have luminosities that are too small. This could be due to a cold 
start (relatively low initial temperatures), a rapid loss of the internal heat, or 
a strong molecular weight gradient that prevents interior regions from cooling 
(Podolak et al. 1991; Hubbard et al. 1995). 

3.7 Do some elements separate from hydrogen at high pressures? 
Where ? 

Helium is strongly suspected of separating from hydrogen in Jupiter and Saturn 
because its lower than protosolar abundance in the atmosphere, and in Saturn 
because without this additional energy source, the planet would evolve to its 
present state in ~ 2 Ga. However, it has not been shown so far that a hydrogen 
helium mixture at Mbar pressures has a critical demixing temperature that is 
above that required in Jupiter and Saturn. 

Helium demixing should occur in the metallic hydrogen region, but it is not 
clear that the critical temperature should decrease with pressure as for fully 
ionized plasmas (Stevenson 1982), or increase with pressure (Pfaffenzeller et al. 
1995). The first scenario would imply the existence of a small inhomogeneous 
region near the molecular/metallic transition as illustrated in fig. El The second 
one would yield a more extended inhomogeneous region. 

Evolution models including the two phase diagrams by Fortney & Hubbard 
(2003) show that in order to reconcile Saturn's age with that of the Solar System 
and the atmospheric helium abundance derived by Conrath & Gautier (2000), 
sufficient energy ~ AMn c gH is required. This implies maximizing H, the dis- 
tance of sedimentation of helium droplets, and hence favors the Pfaffenzeller- 
type phase diagram and the formation of a helium core. 

The question of a phase separation of other elements is still open. It is 
generally regarded as unlikely at least in Jupiter and Saturn because of their 
small abudances relative to hydrogen and the fact that the critical demixing 
temperature depends exponentially on that abundance. 
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3.8 How do the planetary interiors rotate? 

Interior rotation is important because it affects the gravitational moments and 
their interpretation in terms of density profiles (Zharkov & Trubitsyn 1978). It 
is presently not known whether the observed atmospheric zonal flow patterns 
are tied to the planetary interiors or whether they are surface phenomena, with 
the interior rotating close to a solid body with the rate given by the magnetic 
field. Interior rotation affects more significantly gravitational moments of higher 
order. Using extrema set by solid rotation and by a model in which the zonal 
wind pattern is projected into a cylindrical rotation (Hubbard 1982), one can 
show that interior rotation introduces an uncertainty equivalent to the present 
error bar for J4, of the order of the spread in interior models for Jq, and that 
becomes dominant for Js and above. Measurements of high order gravitational 
moments Jg — J14 should tell whether atmospheric zonal flow penetrate into the 
deep interior or whether the deep rotation is mainly solid (Hubbard 1999). 

3.9 What can we tell of the giant planets ' cores ? Are they primor- 
dial? 

Confronted to diagrams such as figs. El and El there is the tendency to think that 
the giant planets cores as well defined, separate entities. It is not necessarily 
the case: first, as shown by fig. 0] solutions with a well-defined central core are 
equivalent to solutions with cores that have been diluted into the central half of 
the planet. Second, convection does not necessarily guarantee the presence of 
globally homogeneous regions, and can efficiently oppose the settling of species, 
as observed in thermohaline convection. Finally, the history of core formation, 
and in particular the epoch at which planctcsimals were accreted and their sizes 
matter (e.g. Stevenson 1985). 

Once formed, the cores of the giant planets are difficult to erode, as this de- 
mands both that heavy elements are (at least partially) soluble in the hydrogen 
helium envelope, and that enough energy is present to overcome the molecular 
weight barrier that is created (Stevenson 1982). However, in the case of Jupiter 
at least, the second condition may not be that difficult to obtain, as only 10% 
of the energy in the first convective cell (in the sense of the mixing length ap- 
proach) needs to be used to dredge up about 20 M e of core material (Guillot 
ct al. 2004). Evaluating whether the first condition is satisfied would require 
knowing the core's composition and its state, but one can nevertherless note 
that the initially high central temperatures (~ 30, 000 K) favor solubility. Such 
an efficient erosion would not occur in Saturn (and much less so in Uranus and 
Neptune) because of its smaller total mass. 

3.10 Do we understand the planets' global compositions? 

This may be the hardest question because it requires tying all the different 
aspects of planet formation to the observations of the atmospheres of the giant 
planets and the constraints on their interior structures. So far, most of the focus 
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has been on explaining the presence of a central core of ~ 10 M® in Jupiter, 
Saturn, Uranus and Neptune. The new interior data suggest that Jupiter's core 
is probably smaller, and that Saturn's may be larger. More importantly, the 
envelopes of all planets appear to be enriched in heavy elements, and this has 
to be explained as well. 

The possibility that Jupiter could have been formed by a direct gravitational 
instability (e.g. Boss 2000) may be appealing in view of its small inferred core. 
However, the enrichment of its envelope in heavy elements is difficult to ex- 
plain within that scenario, given the low accretion rate of a fully-formed Jupiter 
(Guillot & Gladman 2000). 

The leading scenario therefore remains the standard "core accretion" sce- 
nario (Pollack et al. 1996), with the addition that Jupiter, Saturn, Uranus and 
Neptune were closer together (5-20 AU) just after their formation (Levison & 
Morbidelli 2003). While this scenario require core masses ^ 10 M®, the possi- 
bility of an erosion of Jupiter's core is appealing because it would both explain 
the difference in size with Saturn and an enrichment of its envelope. Although 
more limited, a small ~ 2 M® erosion of Saturn's core could provide part of the 
enrichment of the envelope (Guillot et al. 2004). 

The fact that Jupiter's atmosphere is also enriched in noble gases, in par- 
ticular Ar which condenses at very low temperatures (~ 30 K) is still a puzzle. 
Presently invoked explanations include a clathration of noble gases in ices (Gau- 
tier et al. 2001; Hersant et al. 2004), and the delivery of planetesimals formed 
at very low temperatures (Owen et al. 1999). 

Finally, the large enrichments in C and possibly N of the atmospheres of 
Uranus and Neptune probably indicate that a significant mass of planetesimals 
(<; 0.1 M®) impacted the planets after they had captured most of their present 
hydrogen-helium envelopes. Along with the other problems related to this sec- 
tion, this requires quantitative work. 

4 Extrasolar planets 

4.1 Observables 

More than 145 extrasolar planets have been discovered to date (see J. Schnei- 
der's Extrasolar Planets Encyclopedia on http://www.obspm.fr/planets), but 
only those for which a determination of both the planetary mass and radius 
are useful for the purposes of this review. This can only be done for planets 
which transit in front of their star, which, by probabilistic arguments limits us 
to planets that orbit close to their star. I will therefore only be concerned with 
"Pegasi planets", giant planets similar to 51 Peg b and HD209458b (both in the 
constellation Pegasus), with semi-major axes smaller than 0.1 AU. 

Six transiting Pegasi planets have been discovered so far. Their main char- 
acteristics are listed in Table 01 The first one, HD209458b (Charbonneau et al. 
2000; Henry et al. 2000), has been shown to possess sodium in its atmosphere 
(Charbonneau et al. 2002) and to have an extended, evaporating atmosphere 
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Table 3: Systems with transiting Pegasi planets discovered so far 





Age [Ga] 


[Fe/H] 


a [AU] 


T c \ [K] 


M p /Mj 


E p /10 lu cm 


HD209458 a 


4-7 


0.00(2) 


0.0462(20) 


1460(120) 


0.69(2) 


1.02(9) 


OGLE-56 b 


2-4 


0.0(3) 


0.0225(4) 


1990(140) 


1.45(23) 


0.88(11) 


OGLE-113 c 


? 


0.14(14) 


0.0228(6) 


1330(80) 


0.765(25) 


0.77(±|) 


OGLE-132 d 


0- 1.4 


0.43(18) 


0.0307(5) 


2110(150) 


1.19(13) 


0.81(6) 


OGLE-lll° 


? 


0.12(28) 


0.0470(10) 


1040(160) 


0.53(11) 


0.71(12) 


TrES-l f 


? 


0.00(4) 


0.0393(11) 


1180(140) 


0.75(7) 


0.77(4) 



* Equilibrium temperature calculated on the basis of a zero planetary albedo 

a Cody & Sasselov (2002), Brown et al. (2001) 

b Torres et al. (2004), Sasselov (2003), Konacki et al. (2003) 

c Bouchy et al. (2004), Konacki et al. (2004) 

d Moutou et al. (2004) 

°Pont et al. (2004) 

f Laughlin et al. (2004), Sozzetti et al. (2004), Alonso et al. (2004) 

(Vidal-Madjar et al. 2003, 2004). Four others have been discovered by the pho- 
tometric OGLE survey and subsequent radial velocity measurements (Konacki 
et al. 2003, 2004; Bouchy et al. 2004, Pont et al. 2004). One is a result of the 
TrES network survey (Alonso et al. 2004). Present photometric surveys have 
a strong detection bias towards very short periods. Associated to a probability 
of transiting that is inversely proportional to the orbital distance, this shows 
that Table |21 represents only a tiny fraction of planets which may have a low 
probability of existence. 

A crucial parameter for the evolution models is the equilibrium temperature 
T oq = T„^/ R*/2a (assuming a zero albedo, i.e. that all incoming stellar light 
is absorbed by the planetary atmosphere). With values of T cq between ~ 1000 
and 2000 K, the present sample of transiting planets is already quite rich. 

4.2 Observed vs. calculated radii of "Pegasi planets" 

Contrary to the giant planets in our Solar System, Pegasi planets are subject to 
an irradiation from their central star that is so intense that the absorbed stellar 
energy flux is about ~ 10 4 times larger than their intrinsic flux (estimated from 
(jHJl, or calculated directly). The atmosphere is thus prevented from cooling, 
with the consequence that a radiative zone develops and governs the cooling 
and contraction of the interior (Guillot et al. 1996). Typically, for a planet 
like HD209458b, this radiative zone extends to kbar levels, T ~ 4000 K, and is 
located in the outer 5% in radius (0.3% in mass) (Guillot & Showman 2002). 

Problems in the modeling of the evolution of Pegasi planets arise mostly 
because of the uncertain outer boundary condition. The intense stellar flux im- 
plies that the atmospheric temperature profile is extremely dependant upon the 
opacity sources considered. Depending on the chosen composition, the opacity 
data used, the assumed presence of clouds, the geometry considered, resulting 
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temperatures in the deep atmosphere can differ by up to ~ 600 K (Seager & 
Sasselov 1998, 2000; Goukenleuque et al. 2000; Barman et al. 2001; Sudarsky 
et al. 2003; Iro et al. 2004). Because of this problem, and in the framework of 
our simple model, the following discussion will be based on an outer boundary 
condition at lbar and a fixed temperature T\ = 1500 or 2000 K 1 . 

Another related problem is the presence of the radiative zone. Again, the 
composition is unknown and the opacity data are uncertain in this relatively high 
temperature (T ~ 1500 — 3000 K) and high pressure (up to ~ 1 kbar) regime. 
Results from our models are based on opacities from Allard et al. (2001). Other 
calculations using e.g. the widely used Alexander & Ferguson (1994) opacities 
do yield only a slightly faster cooling even though the Rosseland opacities are 
lower by a factor ~ 3 in this regime. 

The resulting mass-radius relations are shown in fig. [23 for T\ = 1500 and 
2000 K, and compared to the observations for the planets listed in Table El For 
each case, an upper limit on the radius is obtained from a pure hydrogen-helium 
composition with Y = 0.25. An ad hoc lower limit comes from a model with a 
15 M ffi central core, and a Y = 0.30 envelope. In both case, the opacity table is 
unchanged. 

Figure ITUI shows that within uncertainties, the measurements for 4 planets 
out of 6 can be explained in the framework of our simple model. However, two 
cases stand out: OGLE-TR-132b appears too small for its age implying that 
it may contain significant amounts of heavy elements in a core or in its deep 
interior. The case of HD209458b is more problematic: the constraints on its age, 
mass, an deep atmospheric temperature that should be ~ 1500 — 2000 K yield 
radii that are about 10 to 20% smaller than measured (Bodenhefmer et al. 2001, 
2003; Guillot & Showman 2002; Baraffe et al. 2003). The fact that the measured 
radius corresponds to a low-pressure (^mbar) level while the calculated radius 
corresponds to a level near 1 bar is not negligible (Burrows et al. 2004) but too 
small to account for the difference. This is problematic because while it is easy 
to invoke the presence of a massive core to explain the small size of a planet, a 
large size such as that of HD209458b may require an additional energy source. 

Bodenheimer et al. (2001) proposed that this large radius may be due to 
a small forced eccentricity (e ~ 0.03) of HD209458b, and subsequent tidal dis- 
sipation in the planet interior. In this case, e > in the energy conservation 
equation ((@J). Because of the relatively limited amount of energy available in 
the (non-circular) orbit and the presumably rapid dissipation (due to a tidal Q 
that is presumably similar to that of Jupiter, i.e. Q ~ 10 5 — 10 6 ), this requires 
the presence of an unseen eccentric companion. The search for this companion 
and a possible non-zero eccentricity of HD209458b is ongoing (Bodenheimer et 
al. 2003). 

A natural possibility may be the stellar flux itself, since transporting to deep 

1 Technically, in order to obtain high entropy initial conditions I use Ti T e q(l-|-L/Leq) 1 ^ 4 , 
but the precise form does not matter as long as L <C £ C q, or equivalently —TidSi/dt -C 
— Ti nt dSi n t/dt where Si n t is the characteristic interior entropy. The 1. equality between T\ and 
Tcq is only a very rough estimate guided by present works on atmospheric models of heavily 
irradiated planets. 
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Figure 10: Mass-radius relation of strongly irradiated planets with ages of 1 Ga 
(upper panels) and 5Ga (lower panels), and 1-bar temperatures equal to 2000 
(left panels) and 1500 K (right panets), respectively. The hashed areas have 
upper and lower envelopes defined by (Y — 0.25, M core = 0) and (Y — 0.30, 
M COIe — 15M(g), respectively. Dotted symbols with error bars indicate known 
objects, plotted as a function of their estimated 1-bar temperatures and ages. 
Planets whose age is uncertain appear in both upper and lower panels. Results 
for non-irradiated planets (dotted lines) are shown for an easier comparison. 
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levels (~ 100 bars or more) only a small fraction of order 0.1% to 1% of the in- 
coming flux would yield a radius that is in agreement with the observations. On 
this basis, Showman & Guillot (2002) proposed that kinetic energy generated 
in the atmosphere due to the strong asymmetry in stellar insolation may be 
transported to deep levels and dissipated there, possibly due to a small asyn- 
chronous rotation and its dissipation by stellar tides. Another possibility evoked 
by the authors that km/s atmospheric winds may maintain the atmosphere into 
a shear-unstable, quasi-adiabatic state, which would force temperatures in ex- 
cess of 3000 K at levels between 10 and a few tens of bars. 

It is puzzling that all other recently announced transiting planets do not 
require an additional energy source to explain their size: this is seen in fig. 1101 
which shows that all planets except HD209458b are consistent with the evolu- 
tionary tracks. 

Is there a consistent scenario explaining all the observations? One possibility 
is that, as proposed by Bodcnheimer et al. (2001), HD209458b indeed has an 
eccentric companion. A second possibility is that their orbital histories have 
been very different. Finally, the planets may well have different compositions. 



4.3 How do tides and orbital evolution affect the contraction and 
cooling of Pegasi planets? 

The small orbital eccentricities of Pegasi planets compared to more distant ex- 
trasolar planets tells us that tides raised by the star on the planet have probably 
played an important role in circularizing their orbits, with a timescale estimated 
at - 1 Ga for a planet at 0.05 AU (Rasio et al. 1996; Marcy et al. 1997). Syn- 
chronisation is expected to occur in only Ma timescales (Guillot et al. 1996), 
maybe much less (Lubow et al. 1997). The tides raised by the planet on the 
star also tend to spin up the star which leads to a decay of the planetary orbit. 
It is interesting to note that, with periods of only ~ lday the three OGLE 
planets lie close to the orbital stability threshold (Rasio et al. 1996), or would 
be predicted to fall into the star in Ga timescales or less (Witte & Savonije 
2002; Patzold & Rauer 2002). 

The energies available from circularisation and synchronisation can be use- 
fully compared to the gravitational energy of the planet (e.g. Bodenheimer et 
al. 2001; Showman & Guillot 2002): 

E circ = * GM * M = 3.6 x 10- m 2 (Hi) (1L) (^Y 1 erg, 
a \0.lJ \MqJ \MjJ VIOR©/ 

E syac = h WA*> = 2.4 x io« (Jf- ) (*L f^rV 
y 2 \0.25J \MjJ Vl0 10 cmy yiO^s- 1 J 

E arav = S = 2.4 x 10 — — r- erg, 

where e is the initial eccentricity, a the planet's orbital distance, M is its mass, 
R is its radius, k is the dimensionless radius of gyration, Aw is the change in 
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the planet's spin before and after synchronisation, and S is approximativcly the 
change in the planet's radius (neglecting any structural changes in the calcula- 
tion of -E S ync and E glav ). E glav is the gravitational energy lost by the planet 
when its radius decreases by a factor ~ 6, or alternatively the minimum energy 
required to expand its radius by the same factor. 

The fact that the three energy sources are comparable imply that very early 
in the evolution, circularisation and synchronisation may have played a role, 
perhaps inducing mass loss (Gu et al. 2004). Once a planet has contracted to 
a degenerate, low 9 state, the gravitational energy becomes large, and circu- 
larisation and synchronisation only have a limited role to play. However, two 
reservoirs can be invoked: the orbital energy of a massive eccentric planet that 
would force a non-zero eccentricity of the inner one (Bodcnhcimer ct al. 2001) 
and the absorbed stellar luminosity in its ability to create kinetic energy in the 
atmosphere (Showman & Guillot 2002). 

A major uncertainty related to these processes and how they affect the plan- 
etary structure is to know how and where energy is dissipated. Lubow et al. 
(1997) proposed that a resonant tidal torque is exerted at the outer boundary 
of the inner convection zone, and that dissipation occurs through the damping 
of gravity waves propagating in the outer stable radiative region. Contrary to 
Jupiter, this may be an efficient process because Pegasi planets have a radiative 
region that extends to great depths. Another possibility is through the excita- 
tion of inertial waves in the convective region, a process that would occur also 
in our giant planets (Ogilvie & Lin 2004). The location of the dissipation is not 
clear, however. If it occurs in the atmosphere, the effect of tides on the evolution 
will be limited, whereas they will have a maximum impact if they occur deep 
into the radiative zone (Guillot & Showman 2002). 

If dissipation cannot reach into the deep interior, the planets will not inflate 
significantly when they migrate to their present location. This would imply that 
HD209458b must have migrated from several Alls to its present location in less 
than ~ 10 Ma (Burrows et al. 2000b). In this framework, one could invoke a 
late migration of the OGLE planets (in particular OGLE-TR-132b) to explain 
their relatively small radius compared to HD209458b. 

4.4 How does the composition affects the structure and evolution? 

It is generally believed that giant planets of the mass of Jupiter should have 
near solar composition and relatively small core masses. However, it may not 
be the case: first, Jupiter is in fact relatively significantly enriched in heavy 
elements. Second, while Jupiter is very efficient at ejecting planetesimals from 
the Solar System, Pegasi planets are unable to do so because the local orbital 
speed (GM./a) 1 ' 2 ~ 150 km s" 1 is much larger than the planet's escape velocity 
(2GM/R) 1 ' 2 ~ SOkms- 1 (Guillot & Gladman 2000). Furthermore, most plan- 
etesimals on low e orbits close to the planet would end up impacting the planet, 
not the star (A. Morbidclli, pers. communication 2004). For this reason, models 
of in situ formation of Pegasi planets generally yield large core masses ~ 40 M ffi 
(Bodcnhcimer ct al. 2000). Pegasi planets should therefore be expected to have 
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very different compositions and core masses, depending on the properties of the 
disk of planctcsimals at their formation, the presence of other planets, and their 
orbital evolution. 
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Figure 11: Evolution of giant planets in terms of radius vs. time, for different 
irradiation levels, and 2 assumed compositions: solar, and 6 times solar. (This 
calculation ignores second order effects as modifications of the adiabatic tem- 
perature gradient and non-linear effects in the opacity calculation, and more 
importantly modifications of atmospheric properties.). 

The presence of a core has a relatively straightforward impact on the evolu- 
tion of giants planets. As shown in fig. El it leads to a much faster contraction 
and a smaller radius at any given age. An enrichment of the envelope both in- 
creases the mean molecular weight and the opacities, with two opposite effects 
in terms of the planet's contraction and cooling. Figure UTI shows that for large 
irradiations (extended radiative zones), the second effect wins and leads to a 
(limited) increase of the planetary radius. However, planets with a larger mean 
molecular weight eventually become smaller. 

The difference in inferred radii between HD209458b and other transiting 
planets could hence indicate that stellar tides play a role in slowing or even 
stalling the contraction of all planets, but that because of different histories, 
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some planets have a large core mass but HD209458b has not. In that framework, 
OGLE-TR-132b would probably need a core of ~ 20 M 9 or more (or the same 
amount of heavy elements in its deep interior) to explain its small radius. The 
large [Fe/H] value measured for its parent star (Table is an indication that 
the planet may indeed have grown a large core. 

4.5 What is the role of the atmosphere for the evolution? 

I have purposedly used a very simple atmospheric model by setting T\ oc T oq = 
cte. Of course, this hides many important complications like opacities, chem- 
istry, gravity dependance, presence of clouds, atmospheric dynamics, depen- 
dance on the incoming stellar flux... etc. These complications partially explain 
differences between several authors (Seager & Sasselov 1998, 2000; Gouken- 
leuque et al. 2000; Barman et al. 2001; Sudarsky et al. 2003; Iro et al. 2004). 
These works yield characteristic temperatures at the base of the atmosphere (i.e. 
where most of the incoming flux has been absorbed) that range from ~ 1700 to 
~ 2300 K. 

However, the largest differences arise from simple geometrical reasons: Be- 
cause these calculations are one-dimensional, some authors choose to model the 
atmosphere at the substellar point, some average the received stellar flux over 
the day-hemisphere (1/2 less flux), and others average it over the entire planet 
(1/4 less flux). This points to real problems: how does the planet reacts to 
this extremely inhomogeneous stellar irradiation, and how do possible inhomo- 
geneities in the atmosphere affect the planetary evolution? 

Without atmospheric dynamics, a synchronous Pegasi planet at ~ 0.05 AU 
of a G-type star would see its substellar point heated to ~ 2500 K or more, and 
its night hemisphere and poles have temperatures ~ 100 K, a clearly unstable 
situation. Assuming synchronisation of the convective interior and a radia- 
tive atmosphere obeying the Richardson shear-instability criterion, Showman 
& Guillot (2002) showed that the atmosphere of Pegasi planets are likely to 
develop km/s winds, but that spatial photospheric temperature variations of 
~ 500 K are likely. Dynamical models using shallow-water equations by Cho et 
al. (2003) also yield latitudinal temperature variations, but predict a surpris- 
ing time-dependent behavior, with a night-side that sometimes becomes hotter 
than the day side. A time-dependant approach of radiative transfer, in which 
the atmosphere is allowed to react to a varying irradiation, shows that a km/s 
rotation indeed yields a ~ 500 K effective temperature variation. It also shows 
that the conditions required for the shallow-water treatment (a relatively long 
radiative timescale) are probably not met in Pegasi planets (Iro et al. 2004). 

As shown by Guillot & Showman (2002), to first order (i.e. neglecting pos- 
sible non-linear behavior due to e.g. opacity temperature dependances and/or 
cloud formation), the cooling with an inhomogeneous boundary condition is 
faster than if the same amount of heat has been homogeneous distributed. This 
is because heat tends to escape more rapidly in regions of low atmospheric 
temperatures. But since the radiative timescale below optical depth unity is 
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approximatively oc P 2 , levels deeper than a few bars tend to homogeneize hori- 
zontally very efficiently, even with a slow circulation (Iro et al. 2004). 

Therefore, there is presently no reason to use for evolution models an at- 
mospheric boundary condition other than that obtained assuming a stellar flux 
averaged over the entire planet. Of course, more work is to be done, as opacity 
variations, the presence of clouds either on the day or night side (depending on 
the kind of circulation) , non-equilibrium chemistry, and possible shear instabil- 
ities and gravity waves damping can all play an important role. 

4.6 Stability and evaporation? 

Because Pegasi planets are so close to their star, the question of their survival 
has been among the first following the discovery of 51 Peg B. Guillot et al. 
(1996) and Lin et al. (1996) independently concluded to a relatively fast con- 
traction of the planet and to its survival based on non-thermal evaporation rates 
extrapolated from Jupiter. These evaporation rates ~ 10 _16 M©a _1 turn out 
to be extremely close to those inferred from observations of HD209458b showing 
the escape of HI (Vidal-Madjar et al. 2003), 01 and CII (Vidal-Madjar et al. 
2004). However, the atmospheric escape problem is more complex than initially 
envisioned, with XUV heating, conduction and gravity waves playing important 
roles (Lammer et al. 2003; Lecavelier des Etangs et al. 2004). 

Generally, a critical question is that of the stability of planets at close orbital 
distances in their young ages (Baraffe et al. 2004; Gu et al. 2004). Figure ITTI 
shows that the cooling timescale is initially relatively long in the case of intense 
irradiation (see also fig. 2 of Guillot et al. 1996) and might lead to a significant 
mass loss in case of a rapid inward migration because of Roche lobe overflow 
(part of the planetary envelope becomes unbound because of the star's gravi- 
tational potential) (Trilling et al. 1998, 2002). Baraffe et al. (2004) find that 
another route may be the strong exospheric evaporation. Below a critical mass, 
the planet would inflate before it can become degenerate enough. However, ei- 
ther the presence of a core and the consequent rapid contraction (see fig. HU|l . or 
an internal cooling associated to the decompression upon mass loss may protect 
the planets from an exponential evaporation. 

5 Conclusion & prospects 

We are just beginning to discover the diversity of giant planets. Already, a 
variety of problems that are particular to one planet or a small ensemble of 
planets have arisen. Given the limited ensemble of objects that we are given to 
study and the rapid evolution of the subject, any attempt to find general rules 
it fraught with risk. Some salient conclusions should however resist the trial of 
time: 

• The giant planets of our Solar System all contain a minimum of 10 M© 
of heavy elements, and even ~ 20 M e for Saturn and probably Jupiter. 
In Jupiter, most of the heavy elements are mixed in the hydrogen-helium 
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envelope. On the contrary, Saturn, Uranus and Neptune appear to be 
significantly differentiated. 



• The envelopes of Jupiter, Saturn, Uranus and Neptune are enriched in 
heavy elements compared to a solar composition, implying that heavy 
elements were delivered cither after the formation (requiring large masses 
in planetesimals because of the low accretion probabilities) or when the 
planets, and in particular Jupiter, were not fully formed. In that case, an 
upward mixing (erosion) of these elements with the envelope is required. 
A third possibility is that these elements were captured in an enriched 
nebula. 

• The demixing of helium in metallic hydrogen has probably begun in Jupiter, 
and has been present in Saturn for 2 — 3 Ga. 

• Like Jupiter and Saturn, the Pegasi planets discovered so far are mostly 
made of hydrogen and helium, but their precise composition depends on 
how tidal effects lead to the dissipation of heat in their interior. 

Improvements on our knowledge of the giant planets requires a variety of 
efforts. Fortunately, nearly all of these are addressed at least partially by ade- 
quate projects in the next few years. The efforts that are necessary thus include 
(but are not limited to): 

• Obtain a better EOS of hydrogen, in particular near the molecular/metallic 
transition. This will be addessed by the construction of powerful lasers 
such as the NIF in the US and the MegaJoule laser in France, and by 
innovative experiments such as shocks on pre-compressed samples. One 
of the challenges is not only obtaining higher pressures, but mostly lower 
temperatures than currently possible with single shocks. The parallel im- 
provement of computing facilities should allow more extended numerical 
experiments. 

• Calculate hydrogen-helium and hydrogen-water phase diagrams. (Other 
phase diagrams are desirable too, but of lesser immediate importance). 
This should be possible with new numerical experiments. 

• Have a better yardstick to measure solar and protosolar compositions. 
This may be addressed by the analysis of the Genesis mission samples, or 
may require another future mission. 

• Improve the values of J4 and Jq for Saturn. This will be done as part 
of the Cassini-Huygens mission. This should lead to better constraints, 
and possibly a determination of whether the interior of Saturn rotates as 
a solid body. 

• Detect new transiting extrasolar planets, and hopefully some that are 
further from their star. The space missions COROT (2006) and Kepler 
(2007) should provide the detection and characterization of many tens, 
possibly hundreds of giant planets. 
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• Improve the measurement of Jupiter's gravity field, and determine the 
abundance of water in the deep atmosphere. This would be possible either 
from an orbiter, or even with a single fly-by (Bolton et al. 2003). 

Clearly, there is a lot of work on the road, but the prospects for a much 
improved knowledge of giant planets and their formation are bright. 
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